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ABSTRACT 



Aims. The characterization of the stellar population toward young high-mass star-forming regions allows to constrain fundamental 
cluster properties like distance and age. These are essential when using high-mass clusters as probes to conduct Galactic studies. 
Methods. NGC 7538 is a star-forming region with an embedded stellar population only unearthed in the near-infrared (NIR). We 
present the first near-infrared spectro-photometric study of the candidate high-mass stellar content in NGC 7538. We obtained H 
and K spectra of 21 sources with both the multi-object and long-slit modes of LIRIS at the WHT, and complement these data with 
sub-arcsecond JHK^ photometry of the region using the imaging mode of the same instrument. 

Results. We find a wide variety of objects within the studied stellar population of NGC 7538. Our results discriminate between a 
stellar population associated to the H ii region, but not contained within its extent, and several pockets of more recent star formation. 
We report the detection of CO bandhead emission toward several sources as well as other features indicative of a young stellar nature. 
We infer a spectro-photometric distance of 2.7±0.5kpc, an age spread in the range 0.5-2.2 Myr and a total mass ~ 1.7x10' Mq for 
the older population. 

Key words. Stars: early-type - Stars: pre-main sequence - Infrared: stars - Galaxy: open clusters and associations: general 



1. Introduction 

Clusters of total mass >10^ Mq are the natural habitat of high- 
mass stars. Whether motivated by the study of the formation, 
evolution of feedback effects of hig h-mass stars themselves 
jBik et all 1201 Ot iMartins et all l2009t iPuga et all l2009b or as 
probes to investigate the physics and chemi stry evolu tion of 
the Galaxy (Messineo et al. 2009; Davies etd] l2009l) . NIR 
spectroscopic studies have been proven powerful to unveil 
and characterise the complete high-mass stellar content within 
massive clusters. 

Although new NIR surveys, such as UKIDSS and VISTA, 
will surely uncover new obscured massive clusters, some 
clusters, akeady known from optical observations, harbour a 
hidden high-mass stellar c omponent only a ccessible at NIR 
wavelengths (e.g. Cyg OB2. lKnodlsededl2000h . 

NGC 7538 (aka Sh 2-158) is a visible Hii region in the 
Perseus spiral arm and a site of active star formation. The early 
detection of several lu minous NIR and f ar-IR sources in the 
vicinity of this region jCampbell & Perss on 1988), hinted to- 
ward a rather massive nature. In particular, NGC 7538 IRS 1 
is a high-mass (~3OM0) protostar with a CO outflow, known 
to power the ultra-compact (UC)H ii region NGC 7538 A. It also 
has an associated linear methanol maser structure, which might 
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trace a Keplerian-rotating circumstellar disk jPestalozzi et al.l 
'2004). Recently, 6.7 GHz methanol masers have been de- 
tected toward the nearby objects NGC 7538 IRS 9 and 
NGC 7 538 S, tracing other y oung and embedded massive pro- 
tostars (" Pestalozzi et al.ll2006l) . 

Although this region has been widely inspected at long wave- 
lengths (mostly in the submillimetre window) a nd even opti- 
cal sp ectroscopy for two stars has been obtained (iRusseil et al.l 
i2007h . only a few detailed near-IR photometric studies of 
the stellar population in NGC 7538 have been conducted 
(Balog et al. 2004; Ojha et al. 2004; McCaughrean et al. 199 1|). 
At subarcsecond resolution, the region breaks down into several 
areas of interest (see Fig. [1): (i) the presumed powering clus- 
ter, centred at the location of IRS 6, (ii) the subcluster sitting 
amidst the bright NIR reflection nebula (IRS 5), (iii) the clus- 
ter located at the southern rim of the optical H ii region (IRS 4), 
(iv) the active region around IRS 1/2, (v) an embedded stellar 
cluster at the location of NGC 7538 S, only detected at wave- 
lengths longward of 3.6 jum, and (vi) the IRS 9 reflection neb- 
ula. Surprisingly, the stellar density detected in the NIR is high- 
est at the souther n rim of the dust bubble that bounds the op- 
tical Hii region dBalog et al.l |2004|) . between regions (iii) and 
(iv). The distance to NGC 7538 has been several times esti- 
mated and values between 2.2 and 3.5 kpc have been reported 
(iMoreno & Chavarria- Kjl986Hlsrael et al.lll973h . 

In this article, we present NIR spectra of a sample of can- 
didate high-mass stars within NGG 7538 together with sub- 
arcsecond NIR photometry. 

We describe the observations and data reduction in Sect. |2] In 
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Fig. 1. Colour composite image of NGC 7538 obtained with LIRIS {blue: J, green: H, red: K^). Overlayed contours depict the 
extension of the Hii region in the DSS/POSSII-F Red map. Stars observed with the LIRIS-MOS mode are labelled in cyan. Object 
#4 was observed with the same instrument in long-slit mode. The green circles and cross are 2MASS point sources of interest. 



Sect. [3] we analyze the spectro-photometric information of a 
sample of candidate high-mass stars. Sect. |4] is dedicated to the 
estimate the distance, age and mass of the cluster powering the 
H II region. Finally, we conclude in Sect.|5]with a brief summary. 

2. Observations and data reduction 

2.1. NIR photometry 

We conducted JHKs imaging observations of the field 
around NGC 7538 with LIRIS (Acosta Pulido et al. 120031: 
iManchado et al.ll2004l) . mounted at the Cassegrain focus of the 
William Herschel Telescope in La Palma. Observations were ob- 
tained on 2006 July 21 as part of the MASGOMAS programme 
that comprised a total of 45 young clusters in the Milky Way (see 
[m arin-Franch et al. (2009) for a description of the programme). 
The average seeing during the imaging observations of this re- 
gion was a.'6Q, a.'63, and a.'64 for J, H, an d K„ respectiv e ly. We 
followed the procedure described in Mar m-Franch et all (l2009l) 
for the reduction and PSF extraction of the point-like sources 
in the field. The saturation limit for the LIRIS photometry cor- 
responds to Ks ~ 10.5 mag. For brighter objects, we comple- 
mented the photometry with the 2MASS point-source catalogue 
values that had their photometric errors determined. 



2.2. NIR spectroscopy 

Follow-up spectroscopic observations of NGC 7538 were ob- 
tained in 2007 September 20, and 2008 August 19 using 
the multi-object spectroscopic (MOS) mode of the LIRIS 
Instrument. The average seeing during both nights was 0'.'9. The 



designed mask comprised 20 slitlets of 0'.'8 in width and 6"in 
length. We used the low resolution grism (A/AA ~ 700) that cov- 
ers the H and K band simultaneously and the medium-resolution 
grisms in the H and K bands with AjAA ~ 2500 and A/AA ~ 1700, 
respectively. The slitlet mask was oriented along a RA.= 142° 
and centred at ajooo^ 23'' 13"' 39', 52000= +61° 28' 60". 
We obtained the observations following a classic ABBA tele- 
scope nodding pattern with an offset of 2'.'5. In the case of two 
slitlets, the spectra appeared contaminated by other sources that 
fell close to one of the nodding positions (position B). Therefore, 
we extracted those specific spectra on only half of the stacked 
frames (at the position A). Observations of the AOV standard 
star HD223386 were obtained in the different configurations in 
order to remove the telluric features due to the Earth's atmo- 
sphere. 

The data were reduced using LIRISDRQa specific data reduction 
package developed under the IRAF environment by the LIRIS 
team. The routines available for the LIRIS MOS mode make use 
of the a-priori mask design information to trace the slitlets and 
their limits in the spatial direction in a uniformly illuminated 
frame. Special care was taken to make the flat-field correction. 
The limited spatial extent of each slitlet may include detector 
defects at comparable or larger spatial scales which are not de- 
tected within individual flat-field slitlets. We circumvented the 
problem by coadding the spectral response information of all 
apertures into a unique monodimensional response, after shift- 
ing by an appropriate offset along the spectral axis. This re- 
sponse was modelled using a high order spline function. The 
pixel-to-pixel correction was determined for each slitlet from the 
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Fig. 2. (Left panel) CMD of the stellar content within the field of view of the left panel. The red squares represent the main 
sequence (MS) of dwarf stars without extinction at a heliocentric distance of 2.7 kpc. The red discontinuous line represents the 
reddening vector for an 09V star, with intervals of Ay = 5 mag indicated by crosses. Dark squares represent objects located at 
a projected distance smaller than 30" from IRS 6. (Right panel) Colour-colour diagram of NGC 7538. Source #36 lies at [4.84, 
2.64], out of the scale of the figure. The dashed lines indicate the direction of extinction according to Rieke & Lebofskv (198 51) with 
overlayed intervals of Ay = 5 magnitudes indicated by crosses. The dash-dotted line represents the T Tauri locus of IMever et al.l 
(Il997h . 



extracted flat-field spectrum divided by the modelled response. 
We performed a standard subtraction of consecutive AB pairs to 
remove the sky background. Individual flatfielding, retracing, re- 
sizing and extraction for each star in every subtracted frame fol- 
lowed. Argon and Xenon arc frames were used to determine the 
wavelength calibration (rms < lA). The resulting spectra were 
coadded weighting by the relative number of counts to provide a 
final spectrum of each star. 

The wavelength calibrated spectra were divided by a normal- 
ized reference spectrum to remove the telluric absorption fea- 
tures. This normalized reference spectrum was generated from 
the observed spectra of the AOV standard star, divided by a Vega 
model spectrum convolved to the correspondent spectral resolu- 
tion. For this, we constructed an IDL routine that iteratively fits 
the components due to the instrument response, an atmosphere's 
transmissivity template and the Vega model spectrum. 



Source #4 was observed with LIRIS long-slit mode in the 
ranges 1.92-2.4yum and 1.52-1.78/im in the nights of 2007 
September 21 and 2008 June 27, respectively. The intermediate- 
resolution grisms provided spectral resolutions of A/AA ~ 1700 
and A/AA ~ 2500. In this case, we traced the spectrum using 
the continuum of the source and applied an optimal weight to 
extract the flux along the spatial direction. The standard star 
HD223386 was observed in the same configurations in order to 
obtain the teUuric correction. 



3. Results 

3. 1 . Colour-magnitude and colour-colour diagrams 

The left and right panels of Fig. |2] show the colour-magnitude di- 
agram (CMD) and the colour-colour diagram (CCD) of the field 
around NGC 7538. In the figure, we have only depicted those 
objects that have LIRIS counterparts in the H and Ks bands and 
in the three channels for the CMD and CCD, respectively. 
The 2MASS source IRS 4 appea rs resolved into tw o components 
in our LIRIS images (see also lOiha et aL I I2OOI . The brighter 
source to the southeast is saturated in our Ks map, while the 
northwest source #2 is one of the brightest in the field. Therefore, 
we use our LIRIS photometry for source #2, while still dis- 
play the 2MASS integrated photometry of these two point-like 
sources in our CMD (IRS 4). IRS 4 is not detected in the 2MASS 
J band survey, therefore, it is important noticing that it shows 
J-H>1.34 and 1.94+0.05, pointing to a strong NIR ex- 

cess. 

This paper aims at identifying the high-mass stellar content of 
a Galactic massive cluster, thus we selected the candidates for 
spectroscopic follow-up using a NIR photometric criterion. We 
looked for reddened bright objects in high stellar density regions. 
The production design specifications of the slitlet mask hinder 
a complete study of extremely crowded regions. A geometrical 
limitation is also imposed by the LIRIS multi-object mode. The 
spectral coverage for each object depends on the position of the 
coiTesponding slit relative to the centre of the chip in the disper- 
sion direction. Therefore, the follow-up sources must be selected 
in a preferred orientation along the sky. 
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Fig. 3. (Left panel) H-band normalised spectra of the 20 sources with MOS (R~2500). A normalised template telluric spectrum is 
shown in the bottom panel. (Right panel) K-band normalised spectra of the 20 sources with MOS (R~1700). A normalised template 
telluric spectrum is shown in the bottom panel. 



Our CCD evidences the strong variations in extinction 
the sources selected for follow-up spectroscopy. 



3.2. Spectral classification 

We present the higher-resolution spectra of the selected high- 
mass star candidates in Fig. [3] and Fig.lD The lower resolution 
spectra are presented in Fig.|5] 

The inspection of the NIR spectra of the 21 objects observed in 
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Fig. 4. (Left panel) H-band normalised spectrum of source #4 obtained with long-slit (R~2500). (Right panel) K-band normalized 
spectrum of the same source with long-slit (R~1700). A normahsed template telluric spectrum is shown below. 



NGC 7538 suggests a classification into the six major groups 
exemplified in Fig.|6l 



Early-type stars: Only 3 objects within our candidate sam- 
ple exhibit atmospheric spectral features in absorption, but none 
of them show traces of Helium. This evidences that sources 
#35, #25 and #17 are B-type stars. Considering the measured 
Br- 11 and Bry EWs, we classify the early-type stars following 
the sch emes described in Hanso n et al. (1996) and .Hanson et ;d] 
(Il998h . The obtained spectral types range between B1.5 and 
B9V. 

IRS 6 and IRS 5 are very bright sources in the field of 
NGC 7538 that we did not include in our NIR s pectroscopic 
follow -up. These objects were already classified by Russeil et al.l 
(l2007h according to the spectral features detected in their blue 
spectra. However, the strength of the Hel and Hell lines on these 
published spectra indicate that, opposite to what the authors 
state in the article, IRS 6 is of earlier spectral type than IRS 5. 
To investigate this, we performed spectroscopic observations in 
September 2009 with ISIS at the William Hersch el Telescope, 
confirm ing that the spectral classification of iRusseil et alj 
(l2007h is swapped for these two particular objects. This revised 
classification solves the disparity in the estimates of the cluster 
distance derived from each star, as it will be later discussed in 
Sect.lO 



0S(1) emission. However, this spectral feature has a nebular ori- 
gin for the earlier case or it is not a dominant feature for the other 
two. 



CO bandhead emission sources: Objects #6 and #14 show 
CO first-overtone emission in their K-band spectra. These 
CO bandheads are emitted by neutral material with tempera- 
tures in the range 2 000-5000 K and densities of ~10"^cm"^ 
(IScoville et alJll97"9h . CO bandhead and Bry e mission are spec - 
tral features commonly displayed by YSOs dBik et alj|2006l) . 
Surprisingly, their //-band spectra show the presence of both 
nebular emission and Brackett features in absorption. These 
two candidate YSOs are fairly isolated: #6 lies in the centroid 
of a trunk-like structure at the rim of the Sh 2-158 Hii re- 
gion in the IRAC maps, whereas #14 belongs to the embedded 
cluster detected at the location of NGC 7538 S or region (v). 
Sources #1 1 and #2, located respectively 13" and 2" away from 
IRS 4, show Bry and CO first-overtone emission in their spectra. 
Contamination by IRS 4, particularly in the spectrum of source 
#2, must be considered. 



Featureless sources: Sources #23, #36, and #18 are identified 
as featureless. The NIR spectra of these objects could resem ble 
those of low-mass Class I objects in (Greene & Ladalll996h or 
even the Herbig Ae stars found bv iRodgersI (120011) . 



H2 emission sources: NIR ro-vibrational emission from H2 is 
detected toward source #63. This source is located in the vicin- 
ity of the active IRS 1-3 region and a midst the strong extended 
H2 em ission reported bv iDavis et al. (1998) and Bloom er et alj 
(1 19981) and attributed to IRS 1/2 ("Kraus et al. 2006). The molec- 
ular hydrogen emission has indeed a thermal origin, pointing to 
the presence of shocks produced by outflows. 
The spectra of sources #60, #1 1 and #2 also show traces of H2 1- 



Brackett+Fe emission sources: Two objects (#26 and #12) in 
the field show a peculiar spectrum with the Brackett series in 
emission and several Fell lines. 

Source #12 shows broad Brl5-10 features (average FWHM of 
220kms-i) and Fell lines at 1.688, 1.742, 2.061 and 2.089jum 
and [Fell] line emission at 1.644 fj.m. For classical Be stars, 
the presence of Fe II emission in the K band, while Mgll 
is not detected, is observed in stars of spectral type B4-B7 
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Fig. 5. H- and K-band normalised spectra of the 20 sources with MOS (R~700). A normalised template telluric spectrum is shown 
in the bottom panel. 



(IClark & Steeld l2000l). A hint of CO bandhead emission is ob- 
served in its K-band spectrum. Due to the presence of forbidden 
line emission, we classify #12 as a Herbig Be. 
The CO bandhead and broa d Bry emission are distinctive fea- 
tures associated with YSOs jBik et al.ll2006h . The spectrum of 
object #26 (aka IRS 3) shows Brl7-10, Fell lines and prominent 
[Fell] emission at 1 .553, 1 .644 and 1 .677 fim. The ratio of [Fell] 



1.553 to 1.644/vm flux for this source (~0.33 ) is indicative o f 
a high-density environment «(, >10 cm ( Ha mann et a 01994 . 
This is not surprising, since #26 appears reddened in our CCD 
by an Ay ~27 mag. Source #26 appears therefore deeply buried 
and located at only 0.2 pc in projecte d distance from the high 
mass protostars IRS 1/2 dBloomer et al. 1998) . The presence of 
[Fell] lines toward #26, together with the detection of H2 emis- 
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Fig. 6. Sample of representative spectra in NGC 7538. B-star: 
B-type star, H2: dominant H2 lines, YSO: Young Stellar Object, 
FL: featureless spectrum, Br-i-Fe: Brackett and Fell emission 
lines, Cool-St: cool star. 



sion in the same spectra, suggest the presence of dissociative J- 
type shocks, tracing flow activity toward #26 and confirming the 
young nature of this source. Source #26 has an associated Hn 
region studied in radio recombina tion lines and detec ted contin- 
uum counterparts at 3.6 and 2 cm jSewilo et al.ll2004) . These au- 
thor's inferred Lyman continuum photon rate corresponds to that 
of an early B-type dwarf. This result is in good agreement with 
our photometry and the dereddened F(HeI 2.058)/F(Bry)=0.073, 
which indi cates temperatures low er than those produced by O- 
type stars jBeniamin et alj|1999l) . We classify #26 as an H 11 re- 
gion powered by a B star with ongoing flow activity. The pres- 
ence of Fell lines may indicate that although much of the emis- 
sion is from an HII region, the object is still rather young and 
surrounded by circumstellar material. 

We note the strong similarity of the H- and K- band spec- 
tra of #26 an d #12 with some LBVs at low excitation (e.g. 
iGeballe et al.l i2000). showing Brackett lines together with per- 
mitted and forbidden Fell lines in emission. However, the low 
luminosity (< IO'^Lq) of these objects, together with the absence 
of Mgll and Sill lines goes against an LBV classification. This 
should be regarded as a warning for missidentifications in LBV 
infrared spectroscopic searches in our Galaxy. 

CO bandhead absorption sources: Seven objects within 
the sample show typical spectral features of cool stars. Our 
class ificatio n of these s ources fol l ows the prescriptions of 
IWinston et alJ (l2009l) and iBik et alJ (I20TO ) for late-type stars. 
This comparative method uses primarily the Mgl, Nal and Cal 
atomic lines as diagnostics for the temperature and the CO first 
overtone absorption bands as luminosity indicator. This clas- 
sification carries an uncertainty of about ±1 spectral subclass. 
Stars #46 and #4 show spectral features compatible with that 
of a KOIII, whereas star #60 is classified as a K2III. Sources 



#55, #39, #58 and #41 have shallower CO absorption bands than 
that of giants and yet deeper than dwarfs, indicating a luminos- 
ity class IV. Therefore, these objects are considered pre-main 
sequence stars (PMS). 

Despite the small offsets between the nodding positions and 
the narrow aperture used for extracting the spectra, nebular H 11 
emission appears strongly in the spectra of sources #58 and #41. 
These objects are located in the line of sight of areas where a 
strong gradient in the nebular emission is detected in our NIR 
maps. Therefore, these lines should be considered with caution. 
A summary of the photometric and spectroscopic characterisa- 
tion of NGC 7538's stellar content is compiled in Table[T] 
The spectral type classification allows the estimate of the extinc- 
tion toward several sources in NGC 7538. The visual extinctions 
that appear in column (8) of Table[T|have been calculated from 
the spectral type classification and corre spondent stellar colours, 
considering the (iRieke & Lebofskvll985h extinction law. The in- 
trinsic colours of O-type dwarfs were taken from lMartins & Ple3 
(l2006l). wh ile th ose for B dwarfs and giants were extracted from 



'Tokunaga' ("2000). In the case of the PMS stars, we considered 
the intrinsic colours reported by Kenvon & Hartmann (1995i). 
Absolute magnitudes are taken from iMartins & Pled ( 120061) in 
the case of the O stars. 



3.3. Infrared excess 

We have complemented our NIR photometry with IRAC/Spitzer 
mid-infrared (MIR) data for this region (Chavarria et al. private 
communication) in order to study the infrared excess toward our 
follow-up sources. Fig. |7] shows the spectral energy distribution 
(SED) of the spectroscopically identified young sources within 
NGC 7538. 

We fitted the LIRIS and IRAC flux distributions with: (i) a 
Kurucz model depicted with a solid line for objects with a 
spectral type classification of according temperature and sur- 
face gravity; (11) a single black-body distribution represented by 
a broken line for those sources whose spectral type could not be 
constrained. Only in the first case, the fluxes were dereddened 
following the Fitzpatrick (1999) parametrization. The total lu- 
minosity for each object was computed integrating the area be- 
low the respective model fit and it is considered a lower limit for 
objects that display an excess. 

Object #63 shows both NIR and MIR excess, although 
the 8.0//m IRAC map appears saturated at this location. The 
large I RAC colour [3.61-[4 .51=0.9 mag hints toward a Class I 
source. iMartins et alJ (l2009l) have shown that the H2-dominated 
sources in the RCW 79 and RCW 120 star-forming regions 
have bright counterparts at 24.0 fj.m. The authors conclude that 
this large infrared excess would be due to the existence of an 
envelope, implying that these H2 sources would be in an early 
evolutionary phase. 

However, the attribute of massive YSO relies in the additional 
presence of a strong IR excess. #6, #14 and #11 exhibit both 
NIR and MIR excesses. Due to its proximity to the bright 
source IRS 4, and saturation in the IRAC maps, we can only 
report the NIR excess of source #2 at sub-arcsecond spatial 
resolution. The photospheric features present in sources #6 
and #14 seem indicative of a more evolved phase, when the 
central star has significantly eroded the disk and photospheric 
features in absorption are present. This idea is supported by 
their IRAC colour-colour classification: Class II in the case of 
#6 and Class I for #14 and #11. However, the fits of the SEDs 
to black-body distributions indicates a decreasing temperature, 
that could be directly related to the different spectral features 
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Fig. 7. SEDs of young sources in NGC 7538. Black dots with error bars are photometric measurements with LIRIS and IRAC. For 
the sources with a spectral type classification and an extinction estimate, the SED has been dereddened. Solid gray lines represent the 
Kurucz model fit of the respective spectral type classification, indicated by the T^ff and lOxlog g labels. The area below the model 
in the range 0.02-1000//m corresponds to the estimated luminosity. Dashed gray lines represent the fit of a black-body distribution 
with the indicated temperature and luminosity. 



observed. 

Objects #23, #36 and #18 posses featureless NIR spectra. #23 
and #18 show both NIR and MIR excess, and a Class II IRAC 
classification, indicating a very young nature dGreene & Ladal 
11997 ). Source #36 is a deeply embedded source (Av~40 mag) 
that displays the NIR colours of a reddened naked photosphere 
without an intrinsic excess. This hypothesis is supported by the 
its non detection in the IRAC maps. We tentatively classify #36 
as an early-type O star whose shallow hydrogen features are not 
detected at the limited SNR of our NIR spectra. 
The SED of #26 (IRS 3) reveals the presence of a strong MIR 
excess that translates into the IRAC colours of a Class I object. 
We have spectroscopically classified source #12 as a Herbig 
Be star. The SED reveals an IR excess that rises up to 8.0yum. 
Its NIR photometry renders a reddening-free parameter 
Q={J-H)-lJQiH-K,) equivalent to -0.29. Herbig Ae/Be 
star candidates cover a range -1.38<Q<-0.22 according to 
[Hernandez et alj (l2005l) . reinforcing our classification. 
Among the identified PMS stars, only source #58 shows a 
significant mid-IR excess. However, the non detection of #55 
and #39 in the IRAC maps at 5.8 and 8.0jum prevents further 
conclusions on the IR excess of these sources. 
Despite classified as OB stars, the SEDs of sources IRS 5 
and #35 exhibit strong MIR excesses (i.e. the typical colours 
of Class I sources). In the case of #35, this effect could be 
explained by the presence of a nearby object (1'.'2 away), that is 
not resolved in the longer wavelength IRAC channels. 



4. Discussion 

4.1. Distance, age and mass 

A n initial estimate o f the distance based on spectro-photometiy 
bv lBlitz etal](ll982h situated NGC 7538 at an hehocentric dis- 
tance of 2.8±0.9kpc. Meanwhile, kinematical esti mates deliv- 
ered a distance for th e same region as far as 3.5 kpc l lsrael et aTl 
Il973h . iRusseil et al] (l26(j7l) determined the spectral type classi- 
fication of two stars in the H ii region: IRS 6 (09V), and IRS 5 
(03V), inferring spectro-photometric distances of 1.6+0.17 and 
4.24+0.29 kpc, respectively. They estimated an average distance 
of 2.27+0. 15kpc for the entire Sh 2-158 region. The most re- 
cent estimate of the distance to NGC 7538 has been obtained 
measuring trigonometric parallaxes of methanol masers, regu- 
larly associated with high-mass st ar forming region s . This study 
yielded a distance of 2.65+{}}^ kpc lMoscadelU et alj (12009 '). 
Our revised optical classification of the O-type sources IRS 6 
and IRS 5 described in Sect. l3.2l renders a new distance estimate 
to NGC 7538. Considering a classification 03V for the IRS 6 
object, we derive a spectro-photometric distance of 2.99+0.5 kpc 
to the cluster powering the H ii region. Likewise, the classifica- 
tion of IRS 5 yields a value of 2.39+0.4kpc. We report a final 
spectro-photometric distance of 2.7+0.5 kpc to this region. 

Despite the large extinction variations in the region that ham- 
pered a robust age determination, Balog et al. (2004) found an 
older generation of stars (~ 4 Myr) and a younger population of 
faint stars closer to ~ 1 Myr. A fraction of 30% of young stars in 
the region around NGC 7538 was reported to exhibit a near-IR 
excess. An upper limit limit to the age of the cluster can be estab- 
li shed from the most massive classified star member. According 
to lSchaerer & de Koteii (Il997h . a 60 Mq star (corresponding to a 
03V star) has a main-sequence lifetime ~2.2 Myr before it starts 



E. Puga et al.: Near-infrared spectroscopy in NGC 7538 



9 







' O.i'MwX ' ft) Myr: 


-2 


L 0.5 Myr \ \ / ; : J 


-1 


1 liHwr — -""""""^ - ' ' \ 

_ / iviyr ^ _ - \ 






\ ) 





















1 


3Myr_jf \ \ 








10Myr x"" ""X 










2 








1 1 1 'V , 1 





4.) 4.0 3.9 3.8 3.7 3.6 
iog (TJ 



Fig. 8. Extinction corrected Hertzsprung-Russell diagram of the 
identified PMS stars within NGC 7538. The open symbols corre- 
spond to PMS stars without an IR excess detected in their LIRIS 
and IRAC photometry, while the filled ones represent those with 
an IR excess detected. The solid lines represent the pre-main- 
sequence isochro nes, while the dotted lines correspond to evolu- 
tionai-y ti-acks bv iDa Rio et al.l (l2009l) . 



the giant phase. Assuming a coeval star-formation event, the old- 
est stellar content in NGC 7538 must be, therefore, younger than 
2.2 Myr 

At the other end of the mass spectrum, we can use the informa- 
tion derived from the PMS stars to constrain the age of the low- 
mass population around NGC 7538. To construct a Hertzsprung- 
Russell diagram (HRD), we use the spectral type to temper- 
ature conversion from Kenyqn&Jiartmann (1995), including 
the overestimates reported bv ICohen & Kuhil 11979) as temper- 
ature error The absolute K magnitudes are then calculated from 
the dereddened apparent brightness and assuming a distance of 
2.7 kpc. The HRD of the PMS stars identified in NGC 7538 is 
shown in Fig. |8] We have depicted in the figure several pre- 
main-sequen c e ev olutionary tracks and isochrones taken from 
lDaRioetal.1 (12009 ). The mass of these PMS stars is estimated 
in the range 2-4 Mq. The comparison with the PMS isochrones 
yields an age range between 0.5 and 2 Myr for the PMS iden- 
tified population. We conclude, thus, that these two estimates 
result in an age range between 0.5 and 2.2 Myr for the powering 
cluster 

Clearly, other massive young stellar populations are also present 
at several locations around NGC 7538. They range from stars for 
which we observe photospheric spectral features, are detectable 
in the optical, but also exhibit very strong IR excess (IRS 5 at 
region (ii)) to objects whose emission is dominated by the pres- 
ence of a circumstellar disk (region (iii)), and finally to nascent 
stars that present on-going outflow activity (IRS 1/2 and IRS 9 
in regions (iv) and (vi), respectively). 

These results support the previous idea of a star-forming 
sequence in the northwest-southeast direction that may be 
smoother than previously considered. The quantitative analysis 
of the individual ages of these populations is beyond the scope 
of this paper. 



Assuming a Salpeter initial mass function (IMF), and extrap- 
olating it down to 0.8 Mq, we can calculate the stellar mass of a 
cluster Yet, this estimate is only meaningful for a coeval stellar 
population. This should be the case for the powering cluster of 
the H II region, with IRS 6 as the ionizing star. Normalising the 
stellar mass distribution by the detection of one 03V star with a 
mass uncertainty between 47 and 64 Mq, we obtain a total mass 
-1.7x10^ M0 for the cluster 

4.2. Spatial distribution 

The second most massive identified star is IRS 5 and, although is 
located amidst the ionised emission, its infrared excess and prox- 
imity to the rim of the molecular cloud suggest that it may not be 
a member of the central cluster, and region (n) in fact comprises a 
younger generation of stars. Among the spectroscopically iden- 
tified B-type stars, only object #25 is located amidst the ionised 
emission, while sources #17 and #35 are farther away. In the 
case of #35, this source is located in the vicinity of the young 
region IRS 1/2 or (iv) and exhibits an infrared excess typical of 
Class O/I objects. 

In fact, our CMD on the left panel of Fig. |2] reveals a deficit 
of early B-type candidates for sources located within a radius 
of 30" from the 03V ionizing star IRS 6 (objects marked by 
squares). To further study this sparseness of early B-type stars 
associated to the cluster, we identify in the CMD B-type candi- 
dates of compatible /riband brightness and H - ^r~0.4. These 
objects appear indicated by diamonds in the left panel of Fig.|2] 
Analogously to our previous analysis of spectroscopically iden- 
tified B-type stars, a few of the B-type candidates show a mid-IR 
excess (Chavarria et al., private communication) suggesting they 
belong to younger pockets of star formation. Other B-type can- 
didates exhibit only photospheric emission and they are located 
to the east and shouth of IRS 6, spreading beyond the boundary 
of the H n region, possibly due to a mass segregation effect. 

We sketch a possible picture in which the cavity created 
by the ionizing stars in the molecular cloud is open in the ob- 
server's direction and part of its stellar population is located in 
the foreground of the neutral molecular cloud, appearing almost 
aligned with other embedded younger st ar-forming regions . This 
hypothesis is supported by the work of 'Balog et al.l (l2004l) who 
found that the stellar density peak of this region at NIR wave- 
lengths is located at the rim of the H ii region. 
This projection effect is particularly important to bear in mind 
when isolating the stellar populations of the different pockets of 
star-formation that may be triggered by the expansion of H n re- 
gions. 



5. Summary and conclusions 

We have reported the spectro-photometric study of a sample of 
candidate high-mass stars in NGC 7538. We have used the LIRIS 
instrument to obtain JHK, photometry and multi-object and 
long-slit spectroscopy in the H and K bands at low and medium 
resolution. Complementary IRAC photometry of the sample has 
aided in the investigation of infrared excesses. 

Besides the optically known O stars IRS 6 and IRS 5, whose 
spectral type classification we have revised, we have pinpointed 
a population of B-type stars located either at or beyond the rim 
of the H 11 region. In view of this, we confirm that IRS 6 is the 
main ionizing source of Sh 2-158. Considering the non detection 
of spectral features and the lack of infrared excess, we suspect 
that object #36 is a deeply embedded early O-type star in the 
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vicinity of IRS 1-3. However, higher signal-to-noise-ratio obser- 
vations are needed to ascertain this hypothesis. 
We have spectroscopically identified four G and K stars still con- 
tracting to the main sequence, and therefore labelled them as 
PMS. These stars are considered to be the precursors to A and F 
stars. 

The observations of IRS 3, aka #26, are consistent with its previ- 
ous classification as a B-type candidate that has already started 
powering a compact H ii region in a high-density environment 
He >10'^cm"^ near the protostars IRS 1-2. The outflow activity 
traced by the [Fell] and H2 molecular emission also support the 
young nature of this object. 

Five of the observed sources have been classified as YSOs, being 
their spectra dominated either by H2 or CO features. Two stars 
display featureless NIR spectra and a prominent infrared excess, 
while object #12 has been identified as a Herbig Be star. In gen- 
eral, all these YSOs are considered to be the precursors of OB 
stars. 

We establish a discrimination between an older population 
of naked stellar photospheres and various regions of more recent 
star formation around NGC 7538. However, the former popula- 
tion overspreads spatially beyond the Hii region likely due to 
a projection effect that difficults the isolation of the population 
associated to the individual pockets of star formation. We find 
that most of the YSOs are concentrated either at the southern 
rim of the H n region or in the area around the infrared sources 
IRS 1-3. This, together with the known on-going formation of 
the protostars in the field, reinforces the idea of a northwest- 
southeast star-forming sequence. However, the study of the de- 
tailed sequence of events and causality between them can not be 
addressed with the presented data. 

We have estimated a distance of 2.7+0.5 kpc, an age range 
0.5-2. 2Myr, and a total mass -1.7x10^ M© for the evolved 
cluster content within NGC 7538. 

Due to the variety of NIR spectral features shown by young 
objects in this field, NGC 7538 is an ideal region to study NIR 
spectral features of YSO for a large range of masses and evolu- 
tionary phases within a determined environment. 
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Table 1. Coordinates, LIRIS photometry and spectral classification of the selected objects within NGC 7538. YSO: young stellar object; FL: featureless spectrum; PMS: 
pre-main-sequence. 



Name 


ID 


R.A. (J2000) 


DEC (J2000) 


J 




H 




Av (mag) 


Class. 


ST 


Notes 


IRS 6 




23 13 34.3 


+61 30 14.3 


8.71 ±0.02'" 


8.21 i 


: 0.04' 


8.00 ± o.or 


5.4±0.4 




03±1V 


09±1'V+ 


IRS 5 




23 13 30.3 


+61 30 11.6 


iO.JO ± 


0.03 


9.82 d 


: 0.04^ 


9.43 ± 0.04"' 


7.6±0.6 




09±1V 


03±1'V+ 




#35 


7S 1 3 47 336 


+61 77 53 74 


18.18 ± 


0.03 


14.62 


± 0.01 


1 7 ■S5 + fl 01 


35.0+0.2 




R1 5-R3V 






#25 


23 1 3 43 820 


+61 29 28 37 


13.41 ± 


0.01 


12.33 


± 0.01 


11 83 + 02 


9 3+0.3 




B2-B7V 






#17 


23 13 31 933 


+61 27 47.40 


13.83 ± 


0.01 


12.38 


± 0.01 


11.59 ± 0.01 


13,6+0.2 




B5-B9V 






#63 


23 13 45.058 


+61 27 53 88 


18.13 ± 


0.05 


15.14 


± 0.03 


13.18 ± 0.01 




YSO 








#6 


23 13 43.498 


+61 29 37.31 


1 1 Q7 -t- 

11.7/ It 


n ni 

U.Ul 


10.87 


± 0.01 


9.73 ± 0.03 




YSO 




neb 




#14 


23 13 44.538 


+61 27 15.19 


14.89 ± 


0.02 


IZ.o 1 


± u.uz 


11.39 ± 0.01 




YSO 




neb Br-abs 




#11 


23 13 30.688 


+61 29 07.47 


17.50 ± 


0.04 


13.78 


± 0.01 


11.15 ± 0.01 




YSO 




H2, neb 


IRS 4 


#2 


23 13 32.219 


+61 29 07.32 


11.40 ± 


0.01 


10.43 


± 0.05 


9.43 ± 0.09 




YSO 




neb 




#23 


23 13 34.240 


+61 29 09.62 


14.13 ± 


0.01 


12.74 


±0.01 


11.78 ±0.02 


_ 


FL 








#36 


23 13 45.380 


+61 28 04.19 


20.09 ± 


0.14 


15.25 


±0.03 


12.61 ±0.01 




FL 








#18 


23 13 57.268 


+61 28 15.63 


15.58 ± 


0.01 


13.15 


±0.01 


11.67 ±0.01 




FL 






IRS 3 


#26 


23 13 43.620 


+61 28 14.00 


16.82 ± 


0.04 


13.67 


±0.03 


11.95 ±0.03 


29.7±0.7 


Hii 








#12 


23 13 33.537 


+61 28 49.80 


15.72 ± 


0.07 


12.92 


± 0.07 


11.10 ±0.04 




Herbig Be 








#55 


23 13 38.488 


+61 28 47.59 


15.98 ± 


0.02 


13.98 


± 0.01 


12.92 ± 0.02 


16.0±0.3 


PMS 


G9±1IV 






#39 


23 13 40.766 


+61 28 28.20 


14.98 ± 


0.01 


13.33 


± 0.01 


12.54 ±0.01 


11.6±0.2 


PMS 


K1±1IV 






#58 


23 13 31.325 


+61 30 22.24 


15.12 ± 


0.01 


13.73 


± 0.01 


12.98 ±0.01 


10.3±0.2 


PMS 


K1±1IV 


neb Bry 




#41 


23 13 33.391 


+61 29 35.77 


14.43 ± 


0.01 


13.15 


± 0.01 


12.56 ±0.01 


7.9±0.3 


PMS 


K1±1IV/III 






#15 


23 13 57.238 


+61 27 47.40 


14.95 ± 


0.01 


12.60 


± 0.01 


11.56 ±0.01 


16.5±0.2 




G8±1III 






#46 


23 13 34.218 


+61 30 17.64 


14.18 ± 


0.07 


13.03 


±0.09 


12.66 ± 0.08 


4.9±1.4 




K0±1III 






#60 


23 13 44.280 


+61 28 10.10 


15.68 ± 


0.01 


13.91 


± 0.01 


13.05 ± 0.03 


11.9±0.6 




K2±1III 


H2 


IRS 7 


#4 


23 13 36.801 


+61 30 39.59 


8.44 ± 0.04' 


7.63 ± 0.04'" 


7.37 ± O.or 


2.6±0.6 




K0±1III 





Photometry fro m the 2MASS Point S ource Catalog. 
Classification bv lRusseil et alj ( |2007|) 
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